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Abstract 

O . This paper presents the results of simulations of the magnetization field ac 

response (at 2 to 12 GHz) of various submicron ferrite particles (cylindrical 
dots). The ferrites in the present simulations have the spinel structure, ex- 
pressed here by Mi_nZn„Fe204 (where M stands for a divalent metal), and 
the parameters chosen were the following: (a) for n = 0: M = { Fe, Mn, 
Co, Ni, Mg, Cu }; (b) for n = 0.1: M = { Fe, Mg } (mixed ferrites). These 
runs represent full 3D micromagnetic (one-particle) ferrite simulations. We 
find evidences of confined spin waves in all simulations, as well as a complex 
Q I behavior nearby the main resonance peak in the case of the M = { Mg, Cu 

^ : } ferrites. A comparison of the n = and n = 0.1 cases for fixed M reveals a 

significant change in the spectra in M = Mg ferrites, but only a minor change 
^ . in the M = Fe case. An additional larger scale simulation of a 3 by 3 particle 

0\ ', array was performed using similar conditions of the Fe304 (magnetite; n = 0, 

M = Fe) one-particle simulation. We find that the main resonance peak of 
cn ' the Fe304 one-particle simulation is disfigured in the corresponding 3 by 3 

l/^ ■ particle simulation, indicating the extent to which dipolar interactions are 

^ . able to affect the main resonance peak in that magnetic compound. 
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1. Introduction 

Ferrites (ferromagnetic oxides) present convenient dielectric and magnetic 
properties for microwave and millimeter-wave applications, considering their 
relatively large magnetic losses and resistivities [l|. It is well known that 
several physical properties at sub-micron scales, such as the size and shape 
of the particles in the system, their composition and concentration, including 
their interactions, are important factors that shape the characteristics of the 
magnetic material in a sensitive manner |2|, y, Ij]. The main interactions 
among these particles are the dipolar (or long range interactions) and the 
spin exchange interactions. The interplay between these interactions often 
lead to novel and complex magnetic phenomena. Therefore, in order to design 
materials appropriate to specific applications, a thorough understanding of 
these phenomena is needed. 

Micromagnetism addresses the study of magnetism at sub-micron scales 
in the continuum approximation, and its main theoretical equation is the 
so-called Landau-Lifshitz-Gilbert equation (LLG) p, la, [Zl, |8[. It describes 
the magnetization vector field dynamics (the local precessional motion of the 
magnetization vector field), including a phenomenological damping term, un- 
der an "effective" magnetic field, representing various interactions amongst 
the spins. Due to the fact that this is highly nonlinear vector partial dif- 
ferential equation, it is generally solved by numerical methods (analytical 
solutions can only be found in very few cases |8|, lol, llOJ). 

Due to the advance of computer capabilities, micromagetic simulations 
have been carried out with increasing validity, elucidating several complex 



magnetic phenomena, but still with many open questions [ll| . In particular. 



studies of the dynamics of confined spin waves in patterned arrays of mag- 



netic particles in thin films |12l . Il3l . Il4j . llSL Il6l . llTj is of great interest and 



is the subject of the present investigation, in which ferrite particles are the 
constituent elements. Although the literature on micromagnetic simulations 



of ferromagnetic or permalloy particles is quite vast (see, e.g., Ref. [16| and 
references therein), possibly due to the fact that such a magnetic material 
is able to support a reasonable range of magnetic structures (specially the 
vortex structure, relevant to magnetic recording systems), the literature spe- 
cific on ferrimagnetic or ferrite particle simulations is still somewhat scarce. 
An inspection of the OOMMF citation list on June 2009 [l8| revealed more 
than 750 papers that have used that simulator, in which only a few of them 



focused on ferrimagnetic particles/films (e.g., Refs. 19|, |20|, l2l|]) or bilayers 



(e.g., Ref. 22|). Zero-field absorption spectra of magnetite cubic particles 
have been reported in Ref. 23|. This motivates our project to systemat- 
ically investigate the magnetization field ac response of various magnetic 
compounds (apart from permalloy) according to several physical properties, 
such as size and shape of the particles, their composition and concentration, 
inter-particle interactions, etc. 

The main purpose of the present work is to study the absorption spectra 
and magnetization dynamics of full 3D micromagnetic simulations represent- 
ing submicron spinel ferrite (e.g., [l|, (2J] and references therein) particles 
(cylindrical dots). Several of these ferrites, already studied for many decades, 
are now being explored in recent advances in nanotechnology, specially in 
spintronics (e.g., Ref. [25J and references therein). It is well known that 
the saturation magnetization of several ferrites can be increased by a proper 
combination with the non-magnetic zinc ferrite. We have attempted to repro- 
duce qualitatively the effect of an addition of zinc content, in order to have a 
picture of its possible contributions to the resulting spectrum. In that case, 
we have focused on a small addition of zinc, which lies in the linearpart of 
the the relation between saturation magnetization and zinc content |l|. This 
first exploration is intended as a basis for a future systematic numerical work 
exploring several material properties of ferrites of various types. 

Another relevant analysis in the present work resulted from the perfor- 
mance of an additional larger scale simulation consisting of a 3 by 3 particle 



array, in a similar fashion to our previous work with permalloy particles [26 
It was performed using analogous conditions of the Fe304 one-particle sim- 
ulation. That larger simulation was performed with the aim of indicating 
the extent to which dipolar interactions are able to affect the spectrum char- 
acteristics in a ferrite patterned film, that is, one formed of closely spaced 
dots. 

This paper is organized as follows. A summary of the time domain mi- 
cromagnetic simulation setups are given in Sec. II. In Sec. Ill, we describe 
the absorption spectra of the simulations and the equilibrium magnetization 
fields and discuss the results, concluding in Sec. IV. 



2. Materials and Methods 

2.1. The Fundamental Equation of Micromagnetism and Spin Wave Phe- 
nomena 

The Landau-Lifshitz-Gilbert equation is a vector partial differential equa- 
tion for the magnetization vector M, defined as the sum of A^ individual 
magnetic moments jlj {j = 1, . . . , N), specified in a elementary volume dV 
at a position vector r within a magnetic system (particle). Being a contin- 
uum limit expression, it assumes that the direction of M varies continuously 
with position |7|. The LLG equation describes the movement of the mag- 
netization field M{f, t) under the action of a external magnetic field [Hext) 
as precession movement of M around an effective magnetic field {Hgff), de- 
fined as Heff = ~f^o^ f."// ■ It is assumed that Eejf embeds the energy 
of an effective magnetic field, which is in turn generally expressed by the 
sum of four fields composing spin interactions of distinct origins, namely: 
Eeff = Eexch + Eanis + -Emag + Ezee (respectively: the exchange energy, the 
anisotropy energy, the magnetostatic or dipolar energy and the Zeeman or ex- 
ternal magnetic field energy). The resulting equilibrium state of that system 
minimizes the total energy. Other physical parameters of the LLG equa- 
tion are: the saturation magnetization, Mg (determined by the temperature, 
here fixed throughout), the gyromagnetic ratio, 7, and a phenomenological 
(Gilbert) damping constant, a. The resulting dynamics is that in which the 
magnetization vector precesses around the -f^e// field, loosing energy accord- 
ing to the damping term, eventually leading to an alignment of M with -ffe//- 
The LLG equation is therefore written as: 



^^^"' ^^ - -7M(r-', t) X #e// - ^M{f, t) X M{r, t) x H^^f . (1) 



dt ' V w ^jj ^_ 



s 



The magnetization dynamics allows uniform and non-uniform (spatially 
varying) precession movement within the system. An oscillating magnetic 
field Hac at a frequency cjq, applied perpendicularly to the magnetization 
field leads to a coupling of M and Hac-, in which the energy will be absorbed 
by the system from the ac field. The ac field couples to uniform (leading 



28, 29 



to main resonance peak) and to nonuniform (spin wave) modes [27 

In the latter case, one notices that exchange and dipolar interactions may 

contribute to the energy of these modes. According to the Kittel's model 



30|, additional resonances will be found at frequencies ujp = Uo + Dk^, where 
D is a. function of the exchange interaction between adjacent spins, and kp is 
the (quantized) wave vector corresponding to a given spin wave excitation. 



It has been noted (e.g., Ref. [31[) that the resonant peaks associated with 
the exchange interactions are found at the left of the uniform resonance 
peak, and that dipolar interactions would be generally independent of the 
size of the system, leading to an interparticle dipolar coupling field (see, e.g.. 



Refs. |12l . |26|). At a reasonably high signal level, they may exist significant 
coupling between the uniform precession mode and spin wave modes, causing 
alterations on the main resonance line (e.g., broadening and lowering). A 
review of confined spin waves can be found in Demokritov et al. [32]. 

2.2. Simulations Setup 

The present work follows the general procedure described in our previous 



work [26|, where a similar analysis of permalloy cylindrical dots has been 
performed, based in the method outhned by Jung et al. [^3]. We have 
used the freely available integrator OOMMF (Object Oriented Micromag- 
netic Framework) |l8|] in order to numerically integrate the LLG equation 
and evolve the magnetization field of the ferrite particles. The particles were 
circular dots with finite thickness, that is, cylindrical dots of 0.5/im of di- 
ameter and 85 nm of thickness. The chosen spinel ferrites for the present 
simulations (Mi_„Zn„Fe204; M, a divalent metal), were the following: (a) 
for n = 0: M = { Fe, Mn, Co, Ni, Mg, Cu }; (b) for n = 0.1: M = { Fe, 
Mg }. Our simulations were performed in a method suitable for a qualitative 
comparison with Brillouin light scattering spectroscopy measurements J33|, 
noting that formally the strength of the Brillouin cross-section differs quan- 
titatively from the amphtude of the absorption spectrum. We describe the 
details of the simulations as follows. 

In order to study the absorption as a function of frequency, an exter- 
nal magnetic field in the plane of the particle was applied, formed by two 
components: a static (dc) magnetic field {Bdc = /io-f^dc) of 100 mT in the y 
direction, and a varying (ac) magnetic field {Bac = /io-f^ac) of small amplitude 
(1 mT) in the x direction: 

S,, = (l-e-^*)5,,,ocos(a;t), (2) 

with the ac field frequency (/ = a;/(27r)) ranging from 2 to 12 GHz, in steps 
of 0.2 GHz (that is, 51 different OOMMF frequency runs were performed for 



each ferrite simulation). We have discretized the time domain of the apphed 
Bac field at intervals of 0.005 ns, which were used as inputs in the "field range" 
record of OOMMF (stepped linearly by the simulator). The simulations 
were run up to 5 ns, resulting in 1000 outputs (dumps) for each of these 51 
frequency runs. An additional simulation involving a 3 by 3 array of particles 
was performed using the same parameters and conditions of the Fe304 one- 
particle ferrite simulation, but running a smaller set of selected frequencies 
around the resonance peak, due to the high computational demand of this 
simulation. 

We list in Table[T]the main global parameters of the OOMMF, which were 
fixed for all sets of simulations. These global parameters were also adopted 
for the larger scale simulation (the 3x3 array). Note that the exchange 
stiffness has a fixed value in all simulations (of the order ~ 10~^^ J/m). In the 
next section we address in more detail this parameter in context of exchange 
length effects. Table [2] lists particular data of the simulations, specifically the 
value of the saturation magnetization and anisotropy constant adopted for 
each ferrite particle. We have extracted the data from Fig. 4.9 of Ref. [l|, 
which presents the experimental values of the saturation magnetization of 
mixed ferrites (in Bohr magnetons) according to the Zinc content {n value). 
Data was also extracted from Refs. [22J and (3J] (c.f. Table A.l). The 
simulations were executed on a 3 GHz Intel Pentium PC running Kurumin 
Linux, taken an average of ~ 28 hours of computation for each set of one- 
particle simulations, whereas the 3x3 array simulation took about two weeks 
to be run. 

2.3. Exchange stiffness considerations 

In numerical micromagnetism, it is important to observe the restriction 
that, in order to obtain accurate results, the value of the computational cell 



size should not exceed the exchange length (see, e.g. Ref. [10|), defined as 
lex = \/2A/(/ioM^), where A is exchange stiffness of the material. Notice the 
stronger dependence of Ig^ to Mg than to A (e.g., a 10% smaller value for A 
with Ms fixed implies a ~ 5% decrease in the resulting l^x, whereas the same 
10% reduction in Mg with A fixed implies a ~ 23% increase in l^x)- The cell 
size here adopted is 5 nm (enough for a meshing of (9(100) magnetization 
cells along the particle's diameter), and therefore materials with Ig^ above 
that limit are in accordance with the present numerical requirements. 



Based on magnetoresistive measurements. Smith et al. [35[ obtained the 



exchange stiffness for Permalloy within 10% error (^NiPe = 1.05 x 10 ^^ J/m). 



This method improves on previous estimates based on spin-wave ferromag- 
netic resonance spectra, which can give discrepancies of a factor of 2 around 
the value A njfs ~ 1 x 10~^^ J/m. From measurements of the domain width, 



Livingston [36| found for Fe-Nd-B magnets Ape-Nd-B = 1.1 x 10~^^ J/m, 
but this method depends on the measurement of the anisotropy constant Kl 
(the quote value for Apc-Nd-B was increased by a ~ 1.5 factor given a new 
measurement of Kl, as mentioned in a note added in proof). Estimates for A 
can also be obtained by a formula that includes the exchange integral J and 
other parameters (see, e.g., ^STJ] for an estimate of the exchange stiffness of a 
nano crystalline NiQ^^Zno,^ ferrite, although uncertainties are not quoted). 

Given the experimental uncertainties, we have decided to adopt an ad- 
equate order of magnitude value for A (such that the resulting l^x is above 
the computational cell of 5 nm), namely A ~ 10~^^ J/m. In particular, we 
have fixed A = 1.2 x 10~^^ J/m, as usually quoted for magnetite (Ref. |l9|. 



but see e.g. Ref. [23| for a quoted value larger by a factor ~ 3). Notice that, 
by the use of a global relation for ferrimagnetic polycrystals |38|, namely, 
A{T) = {kTc/a){l — T/Tc)"^, one can alternatively infer the A{T) value for 
the ferrites from the lattice constant (a) and Curie temperature (T^) with 
< 20% precision. Using this relation and data collected from literature 



rsj 



(see Refs. [39], [40|, [4l|, and |42|), we find that the resulting values of 
A{T) at room temperature for all ferrites in the present work are well within 
~ 20% of the presently adopted value oi A = 1.2 x 10^^^ J/m; in other words, 
reasonably within current experimental uncertainties. 

With the adopted value for A, we see that lex ~ 8.2 nm for the (M = Fe, 
n = 0.1) ferrite (the highest Mg of the set) and lex ~ 33.8 nm for the (M = 
Cu, n = 0) ferrite (the lowest Ms of the set). Hence the latter range for lex 
is above the cell size, in accordance with the numerical requirements. Notice 
that, for the most critical cases (namely, M = Fe ferrites, with the highest 
Ms values) , one could ask how much an error in the corresponding A value is 
allowed for in order to still be in accordance with the numerical requirements 
(considering that the Ms value is correct). It results that a factor of ~ 1/3 
(namely, a decrease in ~ 33% in A) would result in /e^ <~ 5 nm for the M = 
Fe ferrites. We conclude that the adopted value for A is acceptable for the 
present simulations. 

However, it is important to understand how sensitive our simulations are 
to variations in A to the point that the final results could change appreciably. 



We will investigate the effect of lower values of Alj in supplementary simula- 
tions to be discussed in the next section along with the main simulations 

2.4- Calculation of the Spectra 

For the spectra computation, we have also followed the general procedure 



outlined in Ref. |13| (see also |26[). In order to obtain the absorption spectra 
of the ferrite particles, we proceeded as follows, for each simulation. The first 
1 ns of the averaged magnetization vector in the x direction, {M)^{t < 1 ns), 
has been excluded, and the Fourier transform of the remaining time domain 
data, (M)x(l < t < 5 ns), has been calculated. The amplitude of the 
maximum Fourier peak at each frequency run was then selected for providing 
the absorpion at the respective frequency, hence building up the spectrum of 
each ferrite particle. 

3. Results 

3.1. Main Micromagnetic Simulations 

In this section we outline the main results found in the present work. A 
more detailed (qualitative) analysis will be offered in the next section. 

Fig. [1] shows the resulting spline fit absorption spectra of the ferrite one- 
particle simulations for n = 0. It is observed that the ferrites with higher 
saturation magnetization (M^) have their main resonance peaks at higher 
frequencies. It is also possible to notice in each spectrum the presence of 
small amplitude absorption peaks at the left of the main resonance peak; 
these small peaks appear to increase in amplitude for the ferrites with lower 
Mjj. It is already pronounced in the M = Ni case, and results in a "double- 
like" peak in the case of the M = { Mg, Cu } ferrites (which have very similar 
Ms). As already mentioned, resonant peaks associated with the exchange 
interactions are found at the left of the uniform (main) resonance peak (e.g., 
Ref. (3l|)- Hence, the small peaks found in the spectra are probably confined 
spin-wave excitations of the magnetization field of the particles. 

Fig. |2]shows the spectrum of the simulation representing a Fe(i_„)Zn„Fe204 
{n = 0.1) mixed ferrite compared with that oi n = (Fe304). It can be seen 



"'^ Clearly, for a fixed Ms, lower values of A are of more interest than higher values, since 
the latter are "safe" with respect to the numerical requirements for lex, according to our 
considerations. 



that the main peak in the mixed ferrite is shghtly moved to higher frequen- 
cies. Fig. [3] shows the comparative result for the Mg(^_„)Zn„Fe204 mixed 
ferrite. This is a more complex case. Clearly, peak A (n = case), which 
is quite pronounced, lowers in amplitude significantly in the corresponding 
n = 0.1 case (peak C), whereas peak B appears to be modified into peak D, 
which is at a higher amplitude and frequency. 

In Figs, m and E] (top panel), the simulation output "snapshots" of the 
magnetization vector field related to the peaks of interest of the M = {Cu, 
Mg } ferrites are shown. At each peak, the snapshots were chosen (restricted 
to t > 3 ns) at two points of the ac field cycle {ujt = 7r/2 for the snapshot at 
right, and ut = 3tt/2 for the one at left of those figures). The varying pixel 
tonalities of the particle's snapshots correspond to different values of the x 
component of the magnetization field, which was subsampled to show an 
arrow for the average of 9 vectors per cell element. Both simulations show 
similar results due a close Mg value for these ferrites. It is clear that the 
pronounced peak at the left of the spectra (in both these cases) is of a different 
nature from the one at right: in the former peak, the magnetization field in 
the center of the particle is mostly static and aligned with the direction of the 
dc field; the response of the magnetization field is limited to small oscillations 
of the (nonuniform) magnetization near the edges. The corresponding peak 
at right present instead a quasi-uniform behaviour. 

In both panels of Fig. 0, which refer to a comparison between the n = 
and n = 0.1 cases for fixed M = Mg, one is able to contrast the snapshots 
of each of the peaks of interest as a function of the zinc content addition. 
Clearly, the peaks at left (A, C) show a different magnetization field be- 
haviour than the right ones (B, D), as already pointed out in previously. 
Since the zinc addition is here implemented in a qualitative manner, this 
result must be interpreted as a general trend. 

Fig. |6] shows the snapshots of the magnetization vector field at the main 
resonance peak of each of the one-particle ferrite simulations of Fig. [1] (n = 
0). In the case of the M = { Mg,Cu } ferrites, the snapshots were selected 
from the peaks at the right of their spectra (see Figs. H] and for peaks 
B and D), given that they show a similar nature with respect to the main 
resonance peaks of the other ferrites, as already pointed out. It is observed 
a systematic change in the magnetization field response as a function of the 
saturation magnetization of the ferrite (which increases to the left in Fig. [6]). 
This systematic change is revealed in terms of a higher overall amplitude of 
response for higher Mg as well as an increasingly important presence of small 
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oscillations about a nonuniform static magnetization distribution. Fig. [7] is 
a similar figure to the previous one, but presents instead the snapshots of 
the simulations with addition of zinc content (n = 0.1 runs) as compared to 
their n = counterparts. We will discuss these results in more detail in the 
next section. 

Fig. IH] shows the absorption spectrum of the Fe304 3x3 particles (ar- 
ray) simulation superposed to the corresponding single particle simulation. 
Clearly, the main resonance peak of the Fe304 one-particle simulation is dis- 
figured in the corresponding 3 by 3 particle simulation. There are now 3 
resolved peaks approximately within the region of the single main peak of 
the one-particle simulation counterpart, and these peaks increase in ampli- 
tude for higher frequencies, but never reach the same amplitude of the main 
peak of one-particle run. This figure should provide some indication of the 
extent to which dipolar interactions are able to affect the main resonance 
)eak in that magnetic compound. This result is compatible with that of Ref . 



43l | for cylindrical Permalloy 3x3 dot arrays, in which the fundamental 
mode is found to be split into three modes. 

Fig. M shows the corresponding snapshots of the magnetization vector 
field at each of the three peaks identified in the previous figure, concerning 
the 3 by 3 array simulation. Snapshots number 3 (panel at the right of 
that figure) should be compared with that of Fig. [6], M = Fe one-particle 
simulation. There are several issues to be observed in the 3 by 3 array 
simulation, which will be addressed in more detail in the next section. 

3.2. Supplementary Micromagnetic Simulations 



As mentioned on Sec. 12. 3[ we report on additional simulations performed 
in order to evaluate the impact of smaller values of the exchange stiffness 
constant. A, on our results. As explained in that section, it is interesting to 
analyse that impact for the well-known M = Fe ferrite (magnetite). In other 
words, we have artifically lowered the A value for that one-particle ferrite 
model by a factor 1/3 (simulation labeled "SI") and by 20% ("S2"), see Ta- 
ble El Notice that the SI run brings lex ~ 5 nm (cell size), and therefore is 
expected to bring noticeable change in the results (all other parameters re- 
mained fixed). Indeed, as Fig. [TOJ shows, there is a decrease in the amplitude 
of the main peak as A decreases. Otherwise, the resonance frequency and 
other minor modes at the left of the main peak show little variation. This 
suggests that, in addition to the preliminary considerations already expressed 
in Sec. 12. 3^ our results are qualitatively robust. 

10 



In order to verify the sensitiveness of the appearance of the "three peaks" 
found in the 3x3 particles simulation with respect to a change in some specific 
parameter, we have run additional 3x3 simulations with the same parameters 
of the original one, except for a change in some parameter of our choice. Due 
to the fact that these array simulations are computationally demanding, so 
that a fine-grain covering of the parameter space is prohibitive at this time, 
we have limited our analysis to a small set of additional simulations in order 
to infer possible trends. Also, we have limited the simulations to the 5.0 - 
7.0 GHz frequency range, in steps of 0.2 GHz. Table [3] lists the parameter 
changed in these simulations (labelled "S3" to "S7"). 

Fig. [TT] shows the resulting spectra of the additional 3x3 simulations. 
The "Reference" spectrum is that resulted from the original 3x3 run, that is, 
the same as shown in Fig. [8] (the spectrum with "three peaks", as indicated). 
We have also presented a re-analysis of that original simulation by selecting 
the last 2 ns of the remaining time domain data, 3 < t < 5 ns, for the 
calculation of the Fourier transform (instead of the 1 < t < 5 ns data; see 
Sec. 12. 4p . This selects a clearer steady state condition. We see that (top 
panel of Fig. [TT]) the form of the spectrum is practically unchanged, except 
in amplitude, which is decreased. 

A lower value of A ( "S3" ) also produces a smaller amplitude spectrum, 
with the overall form maintained (except perhaps for the first, smaller peak 
at left), which is the same result as the one-particle cases (see Fig. [TD]). A 
higher value for Mg ( "S4" ) results in significant distortion of the reference 
spectrum, namely: a decrease in amplitude of the peaks, and the first, smaller 
peak at left is not seen in the range of frequencies simulated. An anisotropy 
constant Kl set to zero ( "S5" ; Fig. [TTl bottom panel) shows no significant 
change in the spectrum. On the other hand, the spectrum resulting from a 
larger damping constant (factor of 10) misses entirely the three peaks. This 
is interesting in the light of our previous work, where the same damping 
constant was used, and no splitting of the main peak was found for the 
Permalloy 3x3 particles run (see discussions in |26|), although a splitting 



was indicated in Ref. [43||. Our present analysis thus confirm that a larger 



damping parameter possibly explains the difference in the previous results. 



Finally, the spectrum of the "S7" run (model B3 of |26|), where the dots 
"touch" each other, shows a different spectrum as compared to the reference 
simulation (which in turn has 0.122/im of interparticle spacing). However, a 
splitting of the main resonance mode is also visible. 
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4. Discussion 

The elucidation of peculiar features in the absorption spectra of ferrite 
particles must take into account recent interpretations on the nature and role 
of the spin-wave modes. Modes with nodal planes parallel to the magneti- 
zation are associated with high frequency modes, whereas modes with nodal 
planes perpendicular to the magnetization can exhibit frequencies lower and 
higher than the quasi-uniform mode, and their presence mainly depends on 
the number of nodes and the equilibrium between the dipolar and exchange 
interaction effects (see, e.g., Ref. |3l|)- In the present work, we are inter- 
ested in the overall qualitative magnetization field ac response of various 
ferrite cylindrical dots according to the micromagnetic numerical predictions 
for these systems, in order to have a basis for a more detailed subsequent 
investigation. For definiteness, we list here three possible collective ac re- 



sponses of the magnetization field according to the following criteria [12 



• "Quasi-uniform" behavior (QU): the motion of each arrow is approxi- 
mately the same to that of its neighbors, except for the regions around 
the edge of the particle; 

• "Spin- wave" behavior (SW): the arrows exhibit small oscillations about 
a nonuniform static magnetization distribution; 

• "Edge-like" behavior (ED): the magnetization field in the center of the 
particle is mostly static and aligned with the direction of the dc field; 
the response of the magnetization field is limited to small oscillations 
of the (nonuniform) magnetization near the edges - these modes may 
be influenced by the dipolar field coming from another particle placed 
nearby. 

The observed characteristics of the absorption spectra as linked to a visual 
inspection of the magnetization fields at resonant peaks of interest can be 
classified under those criteria, a subject to which we address now. 

A general trend (n = cases) can be seen by comparing the spectra 
of Fig. [1] with the corresponding snapshots of Fig. [61 The nature of the 
resonance peaks based on the appearance of the snapshots can be inferred, 
which we list below: 

• A remarkable feature is that the main resonance peak in all these sim- 
ulations seem to be of a similar nature and follow a systematic pattern, 
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namely: an increase in amplitude response (in the central body of the 
particle) and the presence of small oscillations about a nonuniform 
static magnetization distribution (SW) - both effects as a function of 
a larger saturation magnetization. In the most extreme case, M = Fe, 
one can see clearly the presence of nonuniformity in the magnetic field 
oscillation. 

• The anisotropy constant Ki appears to play a minor (but noticeable) 
role to affect the above-mentioned trend. For example, let us compare 
the M = Co ferrite with its two "neighbours" (in terms of Ms value), 
namely: M = Mn and M = Ni (c.f. Fig. |6]). If one focuses on the pixel 
tonalities (a measure of the magnetization amplitude), the M = Co 
ferrite shows a smaller extent in tonality of the central region (wherein 
the amplitude of the magnetization field is larger) as compared with 
those of the M = Mn (which has a close, but larger value of Ms) and 
M = Ni (lower Ms) ferrites. It would be expected from the above- 
mentioned trend (an increase in amplitude magnetization response in 
the central body of the particle) that such an extent in tonality for the 
M = Co ferrite would be of intermediate size (between the M = Mn and 
M = Ni ones). The fact that this is not observed points to a relatively 
"easier" alignment of the magnetization with the external field in the 
M = Co ferrite case. A more systematic study fixing Mg and varying 
Ki, however, was not performed, and more study is needed to confirm 
these general trends. 

It is specially interesting to observe the snapshots in the cases of the 
"double-peaks" seen in the compounds M = { Cu, Mg } (Figs. H] and O top 
panel), where the pronounced peak at right of the spectrum (in both cases) 
presents an "edge-like" behaviour (ED). Such materials will probably show 
complex spectra in a properly manufactured patterned film, in which "edge- 
like" modes may be significantly influenced by the dipolar field arising from 
another particle placed nearby (depending on the interparticle spacing and 
possibly other factors). We observe that these "edge- like" effects are smaller 
for higher saturation magnetizations. For instance, a comparison between 
the n = and n = 0.1 case (M = Mg; Figs. |5] and [7] ) clearly shows this 
effect (note that the n = 0.1 ferrite has a larger value of Ms than that of 
n = 0, for a fixed M). 

A comparison between Figs. [21 15] and [7] allow us to qualitatively infer how 
peaks of interest possibly morph from one to another as the zinc content 
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is added in the ferrite particle. For the case M = Fe, the nonuniformity 
(SW), as expressed by the pixel tonality distribution, tends to increase in 
amplitude for n = 0.1 (or larger Mg), specially in the central body of the 
particle. Similarly, for the M = Mg case, "edge-like" effects [ED) decrease 
in amplitude and extent for n = 0.1 (compare snapshots related to peaks A 
and C of Fig. [7]). For the quasi- uniform modes (peaks B and D of Fig. [7]) a 
more uniform magnetization field oscillation (QU) in the central the body of 
the particle is found for the n = 0.1 case. These effects qualitatively explain 
the observed transformation of the corresponding spectra. 

A note is necessary at this point. The n = 0.1 simulations were intended 



as preliminar test-cases for a more systematic subsequent work. In Ref. |44 
for example, an experimental study on the overall energy absorption behav- 
ior of Mn-Zn mixed ferrites in the frequency range of 8 to 12 GHz for various 
chemical compositions is presented. A clear understanding of the behaviour 
shown in that work from the point of view of the magnetic absorption dy- 
namics would be desirable. It would be interesting to analyse and compare 
the corresponding numerical predictions with experimental results in order 
to allow for predictions and guidance for specific applicationqj. The present 
test-cases clearly show that such an analysis is feasible. 

The results for the absorption spectrum of the Fe304 3x3 particles (array) 
simulation is very interesting and follows a previous investigation that we 



have performed on similar permalloy arrays [26|. As already mentioned, the 
main resonance peak of the Fe304 one-particle simulation is disfigured in 
the corresponding 3 by 3 particle simulation, resulting in 3 resolved peaks 
approximately within the region of the original single main peak. Such a 
feature is not observed in the permalloy {Ms = 8.0 x 10^ [A/m]) 3 by 3 array 
simulation of our previous work (c.f. Fig. 4, the simulation AO - one-particle 
run - compared with Al - 3 by 3 array run - of that paper, Ref. |26|). 
A reasonable explanation for not finding the splitting of the main mode in 
our previous work is due to a larger damping factor used in that work, as 
suggested here by our supplementary simulations analysis (Sect. 13. 2p . A 



similar splitting was found and discussed in Ref. J43|. There are, however. 



similarities between the present and previous results, which we will address 



^Although the literature on ferrites is extremely vast, it has proven somewhat difficult 
in the course of this work to find adequate papers to which the present results could be 
directly compared. 
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now. 



• The new peaks in the 3 by 3 array simulation never reach the same 
amphtude of the main peak of one-particle run, and this is also true in 
the permalloy simulation of our previous work, although only one peak 
had been observed in that case. As already noted in that paper, the 
decrease in amplitude of the emerging response is probably due to the 
averaging out of the magnetization field over the array particles, which 
show several modes not present in the one-particle case. 

• A very interesting effect (c.f. Fig. |9]), which we attribute to dipolar 
interactions, which was found in our previous work and is confirmed 
here, is the following. Representing the 3 by 3 array as a matrix. A, 
one can observe that the magnetization field of elements in the Aij and 
A3J rows (j = 1,2,3) evolve in symmetric opposition to each other. 
The central row A2J appears to have its magnetization field evolving 
independently of the other two. 

• The central dot in the 3 by 3 simulation has a particularly unique 
behaviour compared to the others in the array. Comparing snapshot 
number 3 (panel at the right of Fig. [9]) with that of Fig. [6l which 
represents the same M = Fe ferrite dot, but completely isolated, we 
notice that the presence of small oscillations about a nonuniform static 
magnetization distribution (SW), seen in the latter simulation, is only 
present in that central dot of the 3 by 3 array. In other words, this 
effect is apparently attenuated in the other dots of the array, which do 
gain instead a more "edge- like" behaviour (ED). 

• Several dots in the snapshots corresponding to peaks 1 and 2 present 
ED behaviour, and this effect appears to be more intense in the dots 
of snapshot number 2 than of number 1, which is consistent with the 
corresponding amplitude of the peaks in the spectrum (c.f. Fig. |H]). 

The effects outlined above must be of a dipolar nature, but the exact 
prediction of the resulting behavior of the magnetization field (anti-) "syn- 
chronism" as a function of the array symmetry or mutual disposition of the 



particles still needs elucidation (see Ref. 26| for several examples and dis- 
cussion) . 
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5. Conclusion 

In the present paper we have reported on a set of 3D niicromagnetic 
simulations of cylindrical dots supposed here to represent subniicron spinel 
ferrite particles excited by an external periodic magnetic field. We have anal- 
ysed the resulting absorption spectra and the magnetization field behavior 
at modes of interest, limited to the timespan covered in the simulations (5 
ns). We have identified the nature of confined spin waves and small oscilla- 
tions of the (nonuniform) magnetization in the absoption spectra through an 
inspection of the magnetization field at extreme amplitudes of the cycle. A 
qualitative analysis of the magnetization field behaviour for all simulations 
was given. 

The absorption spectra of ferrite particles may present complex behavior 
nearby the main resonance peaks, specially in the cases of M = { Mg, Cu } 
ferrites. It is inferred that a significant change in the absorption spectrum 
can be achieved as the zinc content is added in M = Mg ferrites, but this is 
unlikely in the M = Fe case, at least for a change from n = 0ton = 0.1. 
A study of a larger scale simulation of a 3 by 3 particle array with similar 
conditions of the M = Fe one-particle one shows that the resonance peak of 
the one-particle ferrite simulation is replaced by a "triple" peak or otherwise 
disfigured in agreement with Ref. [43[. We confirm our previous result that 
there is indeed a magnetization field (anti-) "synchronism" effect in the array. 
This study permitted us to infer the extent which dipolar interactions are 
able to affect the main resonance peaks in such ferrite particles. 

We aim to perform additional numerical studies to analyse the role of 
confined spin oscillations in various ferrimagnetic particles and arrays with 
different physical conditions in a future work. 
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Table 1: Main parameters set to the OOMMF simulator, fixed for all simulations in the 
present work. 



Simulation Parameter/Option 



Exchange stiffness [J/m] 

Anisotropy Type 

First Anisotropy Direction (x,y,z) 

Second Anisotropy Direction (x,y,z) 

Damping constant 

Gyromagnetic ratio [m/(A.s)] 

Particle thick;ness [nm] 

Particle diameter [/im] ^ 

Cell size [nm] 

Demagnetization algorithm type 

Saturation magnetization [A/m] 



Parameter Value/Option 



f.2 X 10-^^ 
cubic 

(111) 

(10 0) 

0.005 

2.2f X 10^ 

85.0 

0.5 

5.0 

magnetization constant in each cell 

see Table 2 



^ There is a small difference in the case of the 3x3 array simulation. Due to constraints 
in the drawing of the array (bitmap image to be used as input for the simulator), individual 
particles turned out to have 0.552fim of diameter each, with 0.l22fxm of interparticle 
spacing, therefore fitting in an exact square of 1.9fim by 1.9fim. It is necessary that the 
input bitmap size is set as an integer multiple of the cell size. 



f7 



Tabic 2: Particular parameter values for each ferrite simulation 



Simulation (M,n) Ms xlO^ A/m 


Kl xlO^ J/m^; 


Note 


(Fe,0.0) 


5.00 


-1.10 


single particle 


(Fe,0.1) 


5.48 


-1.10 


single particle 


(Fe,0.0) 


5.00 


-1.10 


3 by 3 array of particles 


(Mn,0.0) 


4.14 


-0.28 


single particle 


(Co,0.0) 


3.98 


27.00 


single particle 


(Ni,0.0) 


2.70 


-0.69 


single particle 


(Mg,0.0) 


1.39 


-1.50 


single particle 


(Mg,0.1) 


2.14 


-1.50 


single particle 


(Cu,0.0) 


1.35 


-0.60 


single particle 



Tabic 3: Supplementary simulation data 



Simulation Parameter that has been modified Note 

SI A = 0.40 X 10"^^ [J/m] single particle 

S2 A = 0.96 X 10'^^ [J/m] single particle 

53 A = 0.40 X 10"^^ [J/m] 3 by 3 array of particles 

54 M, 



85 
86 

87 



= u.^u X iu [J/iiiJ o uy o array oi pariicies 

i = 7.5 X 10^ [A/m] 3 by 3 array of particles 

Kl = [J/m^] 3 by 3 array of particles 

damping const. = 0.05 3 by 3 array of particles 

Model B3 of Ref. [26] 3 by 3 array of particles 
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Figure 1: Absorption spectra of single ferrite particle simulations. Curves were arbitrarily 
dislocated for better comparison. Each ferrite spectrum is labelled by its divalent metal M 
(all cases here with n = 0), and organized in order of saturation magnetization, such that 
the upper curve is from the highest Mg . The spectral curves were obtained from spline fits 
of the discrete simulation results (performed at frequencies from 2 to 12 GHz, sampled at 
intervals of 0.2 GHz). Possible spin- wave excitations to the left of main resonance peaks 
can be seen. 
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Figure 2: The spectrum of the simulation representing a Fe(i_„)Zn„Fe204 (n = 0.1) fcrrite 
compared with that of n = (Fe304). 
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Figure 3: Same of previous figure, but for the Mg(]^_„-,Zn„Fe204 fcrrite. Labels A, B, C 
and D mark peaks of interest. 
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Figure 4: "Snapshots" of the magnetization vector field for the CuFe2 04 ferrite, at two 
resonance peaks, as indicated. At each peak, the snapshots were chosen (restricted to 
t > 3 ns) at two points of the ac field cycle (tut = tt/2 for the snapshot at right, and 
wt = 3tt/2 for the one at left). The varying pixel tonalities of the particle's snapshots 
correspond to different values of the x component of the magnetization field (subsampled 
to show an arrow for the average of 9 vectors per cell element). 
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Figure 6: "Snapshots" of the magnetization vector field for the ferrites presented in Fig. 1, 
taken at their main resonance peaks, labehed by its divalent metal M (all cases here with 
n = 0). Satm-ation magnetization increases to the left. In the case of the M = { Mg,Cu } 
ferrites, the snapshots were selected from the peaks at right of their spectra (see Figs. S] 
and[5l and the corresponding explanation in the text). As previously, the snapshots were 
taken at two points of the ac field cycle {ujt = tt/2 for upper snapshot, and ujt = Stt/2 for 
the lower one). 
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of zinc content (n — 0.1 runs) as compared to their n — counterparts. Labels A, B, C 
and D refer to corresponding peaks in Figs. [3] or [5] 
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Figure 8: Comparison between the absorption spectrum of the Fe304 single particle sim- 
ulation (thin line) with that of the 3x3 particles (array) simulation (thick line). Filled 
squares mark the selected frequencies performed for the 3x3 simulation (the associated 
curve is a spline fit to the corresponding data). 
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Figure 9: "Snapshots" of the magnetization vector field for the 3x3 particles simulation 
(n = 0, M = Fe), obtained at the resonance peaks indicated in the previous figure (labelled 
by 1, 2 and 3). The snapshots were taken at two points of the ac field cycle (ujt = 7r/2 
for upper snapshot, and tot = 37r/2 for the lower one). Snapshots number 3 (panel at the 
right) should be compared with that of Fig. El M = Fe one-particle simulation. 
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Figure 10: Comparison of the absorption spectra of a single ferrite particle (M = Fe; 
magnetite), labeled as "Reference" in the figure, with additional simulations "SI" and 
"S2". The latter runs had their stiffness parameter A artifically lowered the by a factor 
1/3 ("SI") and by 20% ("S2") in relation to the reference model. 
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Figure 11: Comparison of the absorption spectra of a 3 x 3 ferrite (M = Fe; magnetite) 
particles simulation, labeled as "Reference" in the figure, with additional simulations "S3" 
to "S7", as listed in Tableland explained in the text. 
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